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We propose a coarse-grained picture to analyze control problems for quantum chaos systems. 
Using optimal control theory, we first show that almost perfect control is achieved for random 
matrix systems and a quantum kicked rotor. Second, under the assumption that the controlled 
dynamics is well described by a Rabi-type oscillaion between unperturbed states, we derive an 
analytic expression for the optimal field. Finally we numerically confirm that the analytic field can 
steer an initial state to a target state in random matrix systems. 
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I. INTRODUCTION 

Controlling quantum systems is one of hot topics in 
physics and chemistry as illustrated in the fields of quan- 
tum information processings P, S Q, ^^'^ laser control 
of atomic and molecular processes [J. As for the lat- 
ter, there have been devised various control schemes: A 
TT-pulse is a simple example to induce a transition be- 
tween two eigenstates 0. As a ^neralization of the tt 
pulse or adiabatic rapid passage we can utilize the 
nonadiabatic transitions induced by laser fields 0]- For 
more than three level systems, STIRAP scheme uses a 
counterintuitive pulse sequence to achieve a perfect pop- 
ulation transfer between two eigenstates Q . When more 
than two electronic states are involved in the controlled 
system, we can use a pulse-timing control (Tannor-Rice) 
scheme to selectively break a chemical bond on a desired 
potential surface by using a pump and dump pulses with 
an appropriate time interval 0. When the controlled 
system has more than two pathways from an initial state 
to a target state, quantum mechanical interference be- 
tween them can be utilized to modify the ratio of prod- 
ucts, which is called coherent control (Shapiro-Brumer) 
scheme [l0| . 

These control schemes are very effective for a certain 
class of processes but are not versatile and ineffective for, 
e.g., multi- level- multi-level transitions we shall consider 
in this paper. There exist several mathematical stud- 
ies which investigate controllability of general quantum- 
mechanical systems 11,12]. The theorem of controllabil- 
ity says that quantum mechanical systems with a discrete 
spectrum under certain conditions have complete control- 
lability in the sense that an initial state can be guided to 
a chosen target state after some time. Although the theo- 
rem guarantees the existence of optimal fields, it does not 
tell us how to construct such a field for a given problem. 

One of the method to practically desi gn an optimal 
field is optimal control theory (OCT) [l^ Il3| | or genetic 



algorithms 0, 0|. We focus on the former in this pa- 
per as a theoretical vehicle. The equations derived from 
OCT are highly nonlinear (and coupled), so we must 
solve them using some iterative procedures. There are 
known some effective algorithms to carry out this pro- 
cedure numerically, however, the field thus obtained is 
so complicated that it is difhcult to analyze the results: 
What kinds of dynamical processes are involved in the 
controlled dynamics? In addition, the cost of the com- 
putation becomes larger if we want to apply OCT to re- 
alistic problems with many degrees of freedom. Several 
efforts have been paid to reduce computational costs; Zhu 
and Rabitz 15j have introduced a non- iterative algorithm 
for the optimal field. 

On the other hand, we know that some chemical re- 
action systems, especially when highly excited, exhibit 
quantum chaotic features |0, i.e., statistical properties 
of eigen-energies and eigen-vectors are very similar to 
those of random matrix systems |l7j . We call such sys- 
tems quantum chaos systems in short. It has been also 
studied how these quantum chaos systems behave un- 
der some external parameters [3 

I23. These sta- 
tistical properties of quantum chaos systems stem from 
multi-level-multi-level interactions of eigenstates, which 
are related to the existence of many avoided crossings 
plf. Hence it is necessary to consider the interaction 
between many eigenstates when we study dynamics in 
such a system. Furthermore, if our purpose is to control 
a Gaussian wavepacket in a quantum chaos system, the 
process also becomes a multi-level-multi-level transition 
because a Gaussian wavepacket in such a system con- 
tains many eigenstates. These are our motivations why 
we treat multi-level-multi-level transitions and want to 
control them. 

This paper is organized as follows. In Sec. ^1 we 
show how quantum chaos systems can be controlled un- 
der the optimal fields obtained by OCT. The examples 
are a random matrix system and a quantum kicked ro- 
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tor. (The former is considered as a strong-chaos-limit 
case and the latter as mixed regular-chaotic cases.) In 
Sec. mil a "coarse-grained" Rabi state is introduced to 
analyze the controlled dynamics in quantum chaos sys- 
tems. We numerically obtain a smooth transition be- 
tween time-dependent states, which justifies the use of 
such a picture. In Sec. IIVI we derive an analytic expres- 
sion for the optimal field under the assumption of the 
CG Rabi state, and numerically show that the field can 
really steer an initial state to a target state in random 
matrix systems. Finally, we summarize the paper and 
discuss further aspects of controlling quantum chaos. 



In actual numerical calculations, we usually solve these 
equations with some iteration procedure [T^ because 
they are nonlinear with respect to \(f>{t)) and The 
optimal field, Eq. lO, is finally given after a local maxi- 
mum of the functional is reached. 

In the following subsections, we numerically demon- 
strate to control multi-level-multi-level transition prob- 
lems in quantum chaos systems: one is a random matrix 
system, and the other is a quantum kicked rotor. 

A. Controlled Random Matrix System 



II. OPTIMAL CONTROL OF QUANTUM 
CHAOS 



We study optimal control problems of quantum chaos 
systems. Our goal of control is to obtain an optimal field 
e{t) which guides a quantum chaos system from an initial 
state \ipi) at t = to a given target state \(ff) at some 
specific time t = T . One such method is optimal control 
theory (OCT), which has been successfully applied to 
atomic and molecular systems Q. 

OCT is usually formulated as a variational problem 
under constraints as follows: We start from the following 
functional used by Zhu-Botina-Rabitz 

J = Jo - a / [e{t)fdt 
Jq 

i-T 



-2Re 



(1) 



The first term in the right-hand side is the squared ab- 
solute value of the final overlap, 



j,^mT)\cpf)\' 



(2) 



The second term is the penalty term with respect to an 
amplitude of the external field e{t). The factor {(p{T)\ipf) 
in the last term is introduced to decouple the conditions 
for the state |0(i)) and the inversely-evolving state jy(i)}, 
which both evolve under the Hamiltonian H[e{t)] pLll3l|. 
The variation of J with respect to \4>{t)) and \x{t)) gives 
Schrodinger's equations, 



(3) 
(4) 



Here we impose the following boundary conditions 

10(0)) = |<^.), \xiT))^\^j). (5) 

Another variation of J with respect to e{t) gives an ex- 
pression for the external field 



e{t) ^Im 
an 



(0Wlx(t))(xWI^|M|0(i)) 



(6) 



The random matrix was first introduced by 
E.P. Wigner as a model to mimick unknown inter- 
actions in nuclei, and has been studied to describe 
statistical natures of spectral fluctuations in quantum 
chaos systems |l7j . Here, we introduce a random matrix 
system driven by a time-dependent external field s{t), 
which is considered as a model of highly excited atoms 
or molecules under an electromagnetic field. We write 
the Hamiltonian 



H[sit)] = Ho+eit)V 



(7) 



where Hq and V are N x N random matrices subject to 
the Gaussian Orthogonal Ensemble (GOE), which repre- 
sent generic quantum systems with time-reversal symme- 
try. The matrix elements of Hq and V are scaled so that 
the nearest-neighbor spacing of eigenvalues of Hq and the 
variance of the off-diagonal elements of V become both 
unity. 

Once we fix the initial state \ipi) and the final state 
Iff), the optimal field e{t) is obtained by some numerical 
procedures for appropriate values of the target time T 
and the penalty factor a. Though there should be many 
situations corresponding to the choice of\(pi) and | </?/), we 
only consider the case where they are Gaussian random 
vectors. It is defined by 



(8) 



where Cj are complex numbers determined from the fol- 
lowing Gaussian distribution. 



P{cj) (X exp I 



(9) 



and 



is an orthonormal basis. |33| We take this state 



because it is typical in a random matrix system. 

We show two numerical examples for a 64 x 64 random 
matrix Hamiltonian: One is the relatively short-time case 
with T = 20 and a ~ 1 shown in Fig. ^ and the other is 
the case with T = 200 and a = 10 shown in Fig. 13 In 
both cases, we obtain the optimal field e(t) after 100 it- 
erations using the Zhu-Botina-Rabitz (ZBR) scheme ^| 
with £{t) =0 as an initial guess of the field. The initial 
and the target state is chosen as Gaussian random vectors 
as mentioned above. The final overlaps are Jq — 0.971 
and 0.982, respectively. 
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FIG. 1: Optimal Control between Gaussian random vectors 
in a 64 X 64 random matrix system by the Zhu-Botina-Rabitz 
scheme with T = 20 and a = 1: (a) the optimal field after 100 
iterations; (b) its power spectrum; (c) the optimal evolution 
of the squared overlap with the target \(4>{t)\ipf)f as well as 
its magnified values near the target time in the inset; (d) 
the convergence behavior of the overlap Jo (solid) and the 
functional J (dashed) versus the number of iteration steps. 



FIG. 2: Optimal Control between Gaussian random vectors 
in a 64 X 64 random matrix system by the Zhu-Botina-Rabitz 
scheme with T = 200 and a — 10: (a) the optimal field 
after 100 iterations; (b) its power spectrum; (c) the optimal 
evolution of the squared overlap with the target, \{(/>{t)\'fif)\^ 
as well as its magnified values near the target time in the 
inset; (d) the convergence behavior of the overlap Jo (solid) 
and the functional J (dashed) versus the number of iteration 
steps. 



One sees that the ZBR scheme is effective enough 
for random matrix systems, i.e., the optimal fields can 
be obtained even for this type of complicated prob- 
lems of multi-level-multi-level transitions. However, it 
seems that the further analysis is difficult because the 
power spectra for the optimal fields, Figs.^b) and|2Ib), 
are very "complex" , i.e., they contain many frequency 
components. |3J| 



B. Controlled Quantum Kicked Rotor 

The kicked rotor (or the standard map) is one of fa- 
mous models in chaotic dynamical systems, and has been 
studied in various situations 0. One feature of its 
chaotic dynamics is the deterministic diffusion along the 
momentum direction. It is also well known that, if we 
quantize this system, this diffusion is suppressed by the 
effects of the wavefunction localization in momentum 
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FIG. 3: Optimal Control in a regular kicked rotor with 
K = 1 and h — 0.3436 by the Zhu-Botina-Rabitz scheme 
with r = 400 and a = 1: (a) the optimal field after 100 it- 
erations; (b) its power spectrum; (c) the optimal evolution of 
the squared overlap with the target \{(f>{t)\ipf)f as well as its 
magnified values near the target time in the inset; (d) the con- 
vergence behavior of the overlap Jo (solid) and the functional 
J (dashed) versus the number of iteration steps. 



FIG. 4: Optimal Control in a chaotic kicked rotor with 
K = 7 and h — 0.3436 by the Zhu-Botina-Rabitz scheme 
with r = 400 and a = 1: (a) the optimal field after 100 it- 
erations; (b) its power spectrum; (c) the optimal evolution of 
the squared overlap with the target, \{<l>(t)\'Pf)f as well as its 
magnified values near the target time in the inset; (d) the con- 
vergence behavior of the overlap Jo (solid) and the functional 
J (dashed) versus the number of iteration steps. 



space p^. 

Here we employ the quantum kicked rotor as a simple 
model of quantum chaos systems. The Hamiltonian of a 
kicked rotor is written as 

i?KR(i) = Y + -COS0 ^{t-nr), (10) 

n— — oo 

where is an angle (mod 27r), p momentum, K a kick 
strength, and r a period between kicks. An external field 



e{t) is applied through the coupling Hamiltonian 

HMt)] = (11) 
where the dipole moment is assumed to be 

fi{e) ^ -cos{e + 660). (12) 

The extra phase 66q is introduced to break symmetry of 
the system. We take SOq = 7r/3 in the numerical calcu- 
lations throughout this paper. The total Hamiltonian is 
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given by 



H[e{t)] ^ HKR{t) + Hi[e{t)]. 



(13) 



For easiness of computation, we impose a periodic 
boundary condition for p as well as 6; the phase space 
of the corresponding classical system becomes a two- 
dimensional torus 22, 23]. In this case, Planck's constant 
is given by ?i = 2ttM/tN, where p = zLMtt defines the 
periodic boundaries in the momentum space, and N is 
the number of discrete points describing 9 and p. In the 
actual calculations, we set t — 1. 

The kicked rotor is often described only at discrete 
time immediately after/before the periodic kicks. In our 
control problem, however, we must represent dynamics 
driven by s{t) between those kicks. Then, we can ap- 
ply Zhu-Botina-Rabitz scheme as usual. According to 
Eq. ijHJl, the optimal external field is given by 

sit) = -±.lm[{mxit)){xitM9)\m] ■ (14) 
an 

Note that, because fj,{9) commutes with the unitary oper- 
ator g-^K cose/h q£ kick, e{t) is obtained as a continuous 
function of time even at the moment of the delta kicks. 

In Figs. 13 and ^ we show numerical results for the 
quantum kicked rotor^^ as in Sec. Ill Al The system 
parameters are chosen to pick up a regular dynamics 
(Fig. |3J) and a chaotic dynamics (Fig. and the oth- 
ers are T — 400 and a = 1. The optimal field after 
100 iterations for the regular case. Fig. |3fa), is rather 
simpler than that for the chaotic case. Fig. EJa). [See 
also Fig. |3|^b) and^b).] This is because more states are 
involved in the latter chaotic process. 

Next we investigate the wavepacket dynamics in phase 
space using the Husimi representation p^ . The initial 
and final states, \ipi) and \(pf), are chosen as minimum 
uncertainty (Gaussian) packets centered at {6i,pi) and 
{9f,pf), respectively. In Fig. Ufa), we show the result 
for the regular case corresponding to Fig. 13 Optimal 
control is achieved for a wavepacket motion within a torus 
with Jo = 0.989. "365 Figure EJb) shows the controlled 
dynamics for the chaotic case corresponding to Fig.^ In 
this case, the wavepacket once spreads all over the phase 
space due to the chaotic nature of the system, but it gets 
together at the target time T with Jq — 0.969. In both 
regular and chaotic cases, the ZBR-OCT scheme works 
well for the quantum kicked rotor j25l |. 



III. COARSE-GRAINED PICTURE 
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FIG. 5: Time evolution of the Husimi distribution for quan- 
tum kicked rotors with h = 0.3436 under an optimal field 
after 100 iterations. The Zhu-Botina-Rabitz scheme was used 
with the penalty factor a = 1 and the target time T = 400. 
From left to right, quantum states immediately after the kick 
at t = 0, 1, 2, 10, 100, 200, 300, 398, 399, and 400 are 
depicted, (a) The parameters are K = 1 (regular case), 
{e„p,) = (1.0,1.0) and (6»/,p/) = (1.0,-1.0); (h) K ^ 7 
(chaotic case), {0i,Pt) = (1.0, 1.0) and (6*/,^/) = (1.0, -10.0). 



where \(f>i) and |02) {Ei and E2) are two eigenstates 
(eigen-energies) of the system, |f2| = |£oMi2|/?i the Rabi 
frequency, ^12 = (0i|/i|02) matrix elements of a dipole 
operator /t, Eq an amplitude of the field, and 9 a certain 
phase parameter. 

In this section, we study the controlled dynamics from 
an initial state \ipi) at t = to a target state \ipf) at 
t = r in a multi-state quantum mechanical system de- 
scribed by Eq. (|7|). By introducing a "coarse-grained" 
picture, which means neglecting highly oscillating terms 
as the case of rotating-wave approximation (RWA) 
and assuming that \ipi) and \ipf) contain many eigen- 
states without any correlation between them, we show 
that the controlled dynamics can be represented as a 
transition between a pair of time-dependent states [2^ . 



If we apply a resonant external field to a two-level sys- 
tem, we can observe a Rabi oscillation. In such a case, 
the quantum state is well described by 

-ie-'^e^''/'^\(f>2)sin\n\t, (15) 



A. Coarse-Grained Rabi State and Frequency 

As shown in Sec. Ill Al the overlap in the controlled 
dynamics rapidly oscillates because the system contains 
many states. To analyze this complicated behavior more 
easily, we introduce the following two time-dependent 
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states, 

\Mt)) = Uo{t,0)\ip,), \xo{t)) = Uoit,T)\ipf) (16) 
where 

«7ofe,ti) = e-'^^°(*^-*^)/'' (17) 

is a "free" -propagator with Hq from t — ti to t2, and T 
is a target time. These states are an analogue of eigen- 
states in the usual Rabi state p5|) , and we try to describe 
the controUed dynamics as a transition from \(f)o{t)) to 
\Xoit)). 



We introduce another quantum state by a linear com- 
bination of the two time-dependent states, 

\m = \Mt))c{t) + |xo(t))5(i) (18) 

where c{t) and s{t) are functions satisfying a normaliza- 
tion condition: 

|c(t)p + |.(i)P = l. (19) 



If we require \4>(t)) to satisfy Schrodinger's equation, we obtain 



ih 



= e{t)V[\MtMt) + \xo{t))s{t)]. 



Multiplying (0o(OI ^'^id (xo(OI from the left gives the following equations for c{t) and s{t) 



in— 



c{t) 
dt \s{t) 



{^o{t)\e{t)V\^o{t)) {MtMt)V\xo{t)) 
iXommiMt)) {xoitMt)V\xoit)) 



c{t) 
s{t) 



(20) 



(21) 



where we have used 



l(0o(Olxo(i))|«l 



(22) 



which is satisfied when \ipi) and \ipf) are random vectors with a large number of elements. 

Our aim is not to solve Eq. 1)21(1 exactly, but to find a coarse-grained (CG) solution by ignoring rapidly oscillating 
terms when the target time T is large enough. If we use the well-optimized field e{t), we expect that the following 
condition 



\{MtMt)v\Mm , \{xo{tMt)v\xom « \{MtMt)v\xom , 



(23) 



are satisfied for T — > oo under the coarse-grained picture. T he validity of this condition will be checked in Sec. 1111 Bl 



Under this condition, we obtain the following simple 
equations 



d_ f c{t)\ _ f tin\ ( c{t) 

dt \s[t) -\hn* j \s{t) 



where 



r.^( {Ut)Ht)v\xo{t)) \^ 



(24) 



(25) 



CG 



is a frequency defined by ignoring rapidly oscillating 
terms. We also expect that Vl has a constant (time- 
independent) value, which will be justified below. Then, 
the boundary conditions c(0) = 1 and s(0) ~ gives a 
solution 

c{t) = cos \n\t, s{t) = -ie"** sin \VL\t (26) 

where e*^ = ri/|f2|. The final expression of the controlled 
dynamics is 

\4>{t)) = |0o(i)> cos \n\t - le-'^lxoit)) sin \n\t. (27) 



Note that this state is interpreted to represent a transi- 
tion between |0o(i)) and |xo(0) O'" that between \(pi) and 
\(pf). Since this is very similar to the usual Rabi state, 
Eq. lO, we caU this state, Eq. "CG Rabi state", 

and the frequency, Eq. "CG Rabi frequency". 



B. Actual Coarse-Graining Procedure 

In the previous subsection, we have introduced the con- 
cept "coarse-graining" (CG) to define the CG Rabi fre- 
quency ri, Eq. (|25|l . In the actual calculations, we carry 
out this procedure by averaging over a certain time in- 
terval. 



(^(i)>CG ^ 



1 



t2~h 



A{t')dt'. 



(28) 



Though this result depends on the choice of ii, t2 in gen- 
eral, we consider that there exists a natural time scale 
where the time averaging is meaningful. In optimal 



7 



control problems, if we choose the target time T large 
enough, we can substitute the range of the integration 
into above expression, i.e., ii = to t2 = T . 

To check when the condition, Eq. H23I) . is fulfilled, and 
when the CG Rabi frequency Q. defined in Eq. (|25|l be- 
comes constant, we introduce the following integrals 

F(t) = f{Mt')W)V\xo{t'))dt' (29) 

g^{t) = [ {c^o{t')\e{t')V\Mt'))dt' (30) 
Jo 

9x{t) = f{xo{t')\e{t')V\xo{t'))dt'. (31) 
Jo 

Though the integrands are rapidly oscillating, a certain 
smoothness can be observed in those integrals, especially 
for F{t). In such a case, we judge that "coarse-graining" 
(CG) is appropriate. Note that F{t) is a linear function 
of t when the CG Rabi frequency is constant. 

FigureElshows |-F(i)|, \g(j,{t)\, and obtained from 

the numerical results in Sec. Ill Al For the case of T = 
20 in Fig. EJa), the values of ^^(i) and Qxit) are small 
compared to F{t), but F{t) cannot be considered as a 
linear function of t. Thus, CG is not appropriate in this 
case. On the other hand, examining the case of T = 200 
in Fig. |HJb), we realize that the condition, Eq. 123|1 . is 
satisfied, and F{t) is regarded as a linear function of t. 
Hence we conclude that CG for random matrix systems 
is appropriate for a rather large target time T, and in 
such a case, the CG Rabi frequency becomes constant. 
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FIG. 6: Absolute values of the functions F{t), g,i,{t), and 
Qxi^) (see the main text) are shown: (a) T = 20 and a = 
1; (b) T = 200 and a = 10. The external fields used in 
these calculations are already shown in Fig.^a) and Fig.|5Ja), 
respectively. 



C. Smooth Transition between Random Vectors 

In Sec. Ill AI we have already obtained the optimal field 
e{t) by the numerical calculation for the random matrix 
systems, Eq. ((TJ. However, only the overlap between 
the time-evolving controlled state \4>{t)) and the target 
state \^Pf) was shown there. In this section, we show the 
overlaps between the time-dependent states defined by 
Eq. (|16() and |0(t)), and find a smooth transition picture. 

In Fig. [7fa), we show the overlap \{4)o{t)\4){t))\'^ and 
|(Xo(i)|'/>(i))P which are obtained from the dynamics 
driven by the same external field shown in Fig. ^a). 
Those curves in the figure are not smooth, and it seems to 
be difficult to approximate them by the CG Rabi state, 
Eq. l|77|l . with a constant fi. In Fig. [Tfb), on the other 
hand, we see a smooth transition from \(l){i{t)) to |xo(0): 
which is induced by the optimal field shown in Fig. I^^a). 
In this case, the dynamics can be well represented by the 
CG Rabi state with a constant Vl. 



IV. ANALYTIC EXPRESSION FOR THE 
OPTIMAL FIELD 

In the previous sections, we have studied the controlled 
dynamics when an optimal field is first given by the ZBR- 



OCT scheme. In this section, in turn, we first assume 
that the dynamics is well approximated by the CG Rabi 
state, and try to derive an analytic optimal field by using 
OCT 2M- 

A. Coarse-Grained Transition Element 

We start from an assumption that optimally controlled 
quantum states are represented by the CG Rabi states, 
i.e., the forwardly evolving state |</'(i)) and the inversely 
evolving state \x{t)) are assumed to be 

\cj){t)) = \(^Q{t))cos\Vl\t~ie~''\xo{t))s\n\n\t (32) 
\x{t)) = -ze^^|(/.o(t))sin|r!|(t-T) 

+ \xo{t))cos\n\{t-T). (33) 

As we have seen numerically in Sec. 1111 CI the optimal 
field induces a smooth transition between |(^o(i)) and 
\Xo{t))- In this section, we employ OCT to study an 
analytic formulation of the optimal field. Substituting 
Eqs. 1(221) and 1(23) into the expression of the optimal 
field, Eq. and after some manipulations, we obtain 

= '^^^^ [e-^' {<^o(t)\V\xm , (34) 
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FIG. 7: The overlaps |(0oW!<^(f)>|^ and \{xo{t)\(l}{t))f are 
shown: (a) T = 20 and a = 1; (b) T = 200 and a = 10. The 
external fields used in these calculations are already shown in 
Fig.^a) and Fig.|5Ja), respectively. 



where \ {(j)Q{t)\xQ{t))\ <C 1 has been used as before. This 
is an analytic expression for the optimal field while the 
value of the CG Rabi frequency fl and the phase param- 
eter 9 have not been determined yet. 

The definition of the CG Rabi frequency, Eq. (|25|1 . is 
used to determine Substituting Eq. and using 
the relation ft = e^^lfll, we obtain 



m = 



sm2\n\T 
4^ 



(35) 



where 



\{Mt)\v\xom' 

+ [e-^'{Mt)\V\xo{t))]' 



CG 



(36) 



is a CG transition element. This equation gives when 
the penalty factor a and the target time T are fixed. 
For a large T, the second term in the right-hand side is 
considered small compared to the first term. In order to 
see this, we represent the initial and final state using the 
eigenstates \(j)k) of as 



(37) 



j k 

with the coefficients Cj and dj. For a large T, we can 



ignore oscillating terms to obtain 

|2 



\{ut)\v\xom 



[{ut)\v\xm? = mT)f 



where 



(39) 



(40) 



becomes small for — > oo when \Lpi) and \^f) are 
random vectors without any special correlation. Thus 
Eq. I|35(l is simplified as 



E 



If the following condition 



> 1, 



(41) 



(42) 



is satisfied, at least one {fl ^ 0) is obtained from 
Eq. I|35() . Using this the final overlap Jq is given by 



Jo = sin^ \n\T, 



(43) 



and the averaged amplitude e of the external field (|34|) is 
calculated as 




mi'dt 



V2h\n\ 

V 



(44) 



In Fig. |S1 we compare the predicted values, Eqs. I|43|) 
and (|44() . with the numerical results for the random ma- 
trix system. Those results agree well each other espe- 
cially for a large T, i.e., the CG picture is valid and useful 
especially for a large target time T.[3^ 



B. Analytic Solution for Perfect Control 

In the ZBR scheme, we must choose a small penalty 
factor a to make the final overlap large enough. In our 
analytical results, if we take the limit a — + 0, we find that 



\n\ 



{2k- l)7r 
2f 



(fc = l,2,...) 



(45) 



satisfies Eq. (|35|l . and then Jo = 1, i-e., a perfect control 
is achieved. Using Eqs. H34|) and H35|l . the optimal field 
for the perfect control in the small a limit is obtained as 

<t) = ^^V4^"^ ^^ [e'^'{Mt)\V\xom (46) 

where 6 can be determined by a normalization condition 
as 



:^0 ^ {MT)\'Pf) 



(47) 



9 



(a) 




1.0 



10 100 
Target Time 



(b) 



0.05 - 
0.04 - 
0.03 - 
0.02 - 
0.01 

0.00 ^ 




1 



10 100 
Target Time 



FIG. 8: (a) The final overlap Jo = \{(t>(T)\ip f)\^ and (b) 
the averaged field amplitude e for a 64 x 64 random matrix 
system are shown as a function of the target time T. Marks 
( X ) represent the numerical results by the Zhu-Botina-Rabitz 
scheme. Solid curves represent our analytic results under the 
assumption of the CG Rabi state. 



O 0.4 - 




FIG. 9: The coarse-grained Rabi oscillation induced by 
the analytical external field for perfect control is shown for 
the case fc = 3 in Eq. 14611 . The solid curve represents 
\{Xo{t)\(j){t))\^ , and the dashed |{</>o(t)|</>(t)>p. The initial and 
the target states are Gaussian random vectors in a 256 x 256 
GOE random matrix Hamiltonian system. 



This field is expected to be the optimal field which steers 
the quantum state \ipi) at i = to at t = T, as well 
as it induces a CG Rabi oscillation between |(/>o(i)) ^nd 
IXo(O)- Note that the penalty factor a does not appear 
in Eq. (|46|l . so this is different from other non- iterative 
optimal fields discussed in |0| . 

We next examine when and how the analytic opti- 
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1000 



FIG. 10: The target-time dependence of the final overlap Jo 
by the analytic optimal field with fc = 1 is shown. The residual 
probability 1 — Jo from the perfect control Jo = 1 is depicted 
for various matrix sizes A'^ of GOE random Hamiltonians. The 
initial and the final states are Gaussian random vectors. 



mal field works for a random matrix system (256 x 256 
GOE random matrix). Figure|51demonstrates the coarse- 
grained Rabi oscillation induced by the analytic field, 
Eq. (|46|l . with fc = 3, where smooth oscillations of 
\{(l)o{t)\(l){t))\'^ and \{xa{t)\(t>{t))?' are observed. The ini- 
tial and the target states are both Gaussian random vec- 
tors with 256 elements. This result shows that the field 
actually produces the CG Rabi oscillation in the random 
matrix system. 

Finally, in Fig. llUI we show the performance of the an- 
alytic field, Eq. , for the same type of control problem 
with various matrix sizes. The abscissa and the ordinate 
are the target time T and the residual probability 1 — Jqj 
respectively. This result shows that the final overlap Jo 
approaches unity, i.e., the perfect control is achieved as 
the target time and the matrix size become both large. 



V. SUMMARY AND DISCUSSION 

We have studied optimal control of random matrix 
systems and a quantum kicked rotor as examples of 
quantum chaos systems. Using the ZBR-OCT scheme, 
we numerically achieved almost perfect control for the 
above systems where the initial state and the 

target state \ipf) are random vectors (except the case 
of a quantum kicked rotor with K = 1). However, the 
optimal fields and the overlap |(0(t)|(/3/)p thus obtained 
are too complicated to be analyzed as shown in Figs. ^ 
121 El and 01 On the other hand, as shown in Fig. [3 
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the overlaps | (0o (01 P and |(xo(*)|</'(0)P a-^'e rather 
smooth where |0o(O) (Ixo(i))) represents a free forward 
(backward) evolution of the system, so we can introduce 
coarse grained concepts: a CG Rabi state and a CG Rabi 
frequency. The CG Rabi state is an analogue of a usual 
Rabi state but it describes a transition between |0o(O) 
and Ixo(i)) as in Eq. (j^Zj). The CG Rabi frequency 
is defined by ignoring rapidly oscillating terms as in 
Eq. if^ . We applied this picture to OCT and obtained 
an analytic expression for the optimal field, Eq. H46|l . 
We also numerically confirmed that the analytic field 
actually works in controlling random vectors when the 
target time and the matrix size are both large enough. 

In closing, we discuss future directions of this study: 
(a) We mainly studied strong-chaos limit cases as de- 
scribed by random matrix Hamiltonians, and applied the 
coarse grained ideas to them. Thus, the next problem 
should be addressed on less chaotic cases as described 
by banded random matrix Hamiltonians. A quantum 
kicked rotor with a small K will be a good example 
for that purpose |23|. (b) The other interesting prob- 



lem is the semiclassical limit of the controlled dynamics. 
Though we have shown that quantum chaos systems can 
be controlled, we don't know its semiclassical behavior 
since there are many difficulties in taking the semiclas- 
sical limit ?i — > 0. There are, on the other hand, many 
works studying chaos control in classical mechanics, and 
there are some examples utilizing stochastic features of 
phase space in "targeting" problems j22j . In this re- 
spect, it is strongly desirable to study chaos control from 
semiclassical points of view 

laSlp. (c) In connection 

with quantum information processings, control of quan- 
tum entanglement in quantum chaos systems [s^ will be 
another interesting subject to be pursued. 
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